arXiv: 1501.03305vl [physics.bio-ph] 14 Jan 2015 


Autonomous movement of chemically powered vesicle 

Shivam Gupta 1 , Sreeja K. K . 2 and Snigdha Thakur 1 
1.Department of Physics, Indian Institute of Science Education and Research Bhopal, Bhopal, India 
2. Department of Physics, Indian Institute of Technology Madras, Chennai, India 

A mechanism for self propulsion of deformable vesicle has been proposed, vesicle moves by sensing 
the self-generated chemical gradient. Like many molecular motors they suffer strong perturbations 
from the environment in which they move as a result of thermal fluctuations and do not rely on 
inertia for their propulsion. Motion of the vesicle is driven by an asymmetric distribution of reaction 
products. The propulsive velocity of the device is calculated as well as the scale of the velocity 
fluctuations. We present the simulation results on velocity of vesicle, reorientation of vesicle, and 
shape transformation of vesicles. 


I. INTRODUCTION 


Nondiffusive transport of objects such as vesicles, col¬ 
loids and aggregates plays a fundamental role in biolog¬ 
ical systems. What is most important for these small 
length scale objects for motility is some degree of asym¬ 
metry. Several mechanisms that can lead to such mo¬ 
tion with the help of asymmetry have been predictecP®. 
One important mechanism responsible for the propulsion 
is the conversion of chemical energy into directed mo¬ 
tion. At cellular level, nano-sized biological molecules 
are transported with remarkable efficiency by biologi¬ 
cal motors such as kinesin through chemically induced 
symmetry-breaking conformational changesIn this 
paper we report the autonomous propulsion of vesicle 
using the above, as an underlying mechanism. 


Inspired by nature synthetic structures have been de¬ 
signed that uses chemical, light or other form of en¬ 
ergy to perform directed motion . 4 Two different strat¬ 
egy are followed for the construction of such motors. In 
one biology is copied and artificial molecular machines 
are constructed such that they use chemical energy to 
drive the conformational changes and hence perform di¬ 
rected motion 8 9 . For the other route the propulsion rely 
on asymmetric chemical reactivity instead of asymmetric 
conformational dynamics. The most widely studied ex¬ 
ample of this type of motors are bimetallic nanorods 10 11 , 
Janus particled^, sphere dimer motors 13 ^ polymers 16 
The above cited example of physically asymmetric self- 
proplled objects assume that the particle shape is un¬ 
changed during the motion. However, in reality many 
self-propelled objects may change their shape depending 
on the velocity, environment or interaction with other ob¬ 
jects. Recently there have been reports of autonomous 
propulsion for deformable objects with propulsive ef¬ 
fect created by the physical asymmetry ? 17 22 Ohta and 
coworkers have theoretically investigated the individual 
and collective moti on of deformable objects exhibiting 
self-propulsion . 19 23 They have found a bifurcation from 
a straight line to a circular motion in case of an isolated 
particle. In a different work Giardini and coworkers have 
shown the propulsion of large artificial lipid vesicles due 
to the forces generated by actin in a motile comet tail . 22 
Another very recent work by Wilson and group have ex¬ 


perimentally demonstrated the autonomous movement of 
platinum-loaded stomatocytesP^ These polymeric ana¬ 
logue of liposomes are created by the diblock copolymers 
and exhibit remarkably high average directed velocity 
even in presence of small amount of fuel- which is why 
they can be seen as ideal drug delivery devices. 


Here, we illustrate the self-propulsion of a chemo- 
mechanically active vesicle. A molecular description of 
chemically powered directed motion was obtained in the 
triangulated vesicle model with localized asymmetric cat¬ 
alytic reactions. Vesicle consists of linked catalytic (C) 
and non-catalytic (N) vertices, which then is immersed 
in a solvent containing reactive A species. An irreversible 
reaction A + C -A B + C takes place at all the catalytic 
vertex, which in turn gives rise to nonequilibrium gra¬ 
dient of B species. The combination of nonequilibrium 
gradients and potential asymmetry gives rise to the di¬ 
rected motion, similar to nanodimersP 4 Here we discuss 
the factors that control the vesicle velocity and show how 
vesicles can be designed to have velocities that are larger 
than the thermal fluctuations, leading to easily observed 
directed motion. Such directed motion is essential for the 
applications where such propelled vesicles may possible 
be used as a carrier of target drug delivery or to perform 
specific task. 


The paper is organized as follows. In Section [TT] we de¬ 
scribe the particle-based mesoscopic model for the vesicle 
motors and solvent. Section [HI] presents various vesicles 
shape obtained in our model and their propulsion be¬ 
havior. The efficiency of the vesicle motors of different 
shapes has been discussed in Section [IV| Finally, the con¬ 
clusions of our study are presented in Section |V| 


II. SIMULATION MODEL 


Self-propulsion of the fluid vesicle is simulated using 
a mesoscopic approach, which combines a particle based 
hydrodynamics model for solvent^ 2 ^ and a mesoscald^, 
dynamically triangulated surface model for the mem- 
brand^P Similar model has been used previously to 
stud y the dynamics of fluid vesicle in capillary and shear 
flows.™ 
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Triangulated-surface model for the vesicle 

The vesicular membrane is modelled using dynam¬ 
ics triangulation method in which the continuum sur¬ 
face is replaced by a descretized surface defined by ver¬ 
tices, triangles and links. A spherical topology with the 
triangulation thus consists of N v vertices connected by 
Nl = 3(N V — 2) links, Nt = 2{N V — 2) triangles. Such 
description of the membrane is extremely useful when 
only the gross features of the structure and interactions 
are important. The curvature elastic energy of the vesi¬ 
cle responsible for the shape and fluctuations is given by 
Canham Helfrich ModeP? 

H c = | J(c i + c 2 ) 2 da (1) 

where, c\ and c 2 denote the two principle curvatures and 
k is the bending stiffness of the membrane. For the tri¬ 
angulated surface we use a discretized form of this Hel¬ 
frich Hamiltonian proposed by Itzykson 29 31 . 

Unlike the Monte Carlo simulations of triangulated 
membrane we use Stillinger Weber type potential® for 
bonding ([/&) the vertices and also for excluded volume 
interaction (U ex ). The attractive bonding potential be¬ 
tween the connected vertices is give by 

f ^xpy/^o—«•)] ifr>rc0 

Mr)= 0 if r < r c0 ( 2 ) 


where r = |ry| is the bond length between vertices i and 
j. The repulsive excluded volume interaction is incorpo¬ 
rated by: 


U ex (r) 


7 exp[l/(r-r cl )] jf r < 

T T min 

0 if r > r c \. 


(3) 


We use Vmin = 0.67ao and r max = 1.33ao as the mini¬ 
mum and maximum possible distance between vertices. 
The time evolution of the system is governed by hybrid 
molecular dynamics-multiparticle collision (MD-MPC) 
dynamical 

The triangulated membrane acquires its lateral fluidity 
from the bond flip mechanism, where the tethers can be 
flipped between the two possible diagonals of two adja¬ 
cent triangles at time interval tbf - 3 ^ Flipping between 
catalytic and non-catalytic part of the vesicle is not al¬ 
lowed. 


Multiparticle collision dynamics of the solvent 

To bridge the large length and timescale gaps in the 
vesicle dynamics a mesoscopic model for the solvent is re¬ 
quired. In multiparticle collision dynamics, the solvent is 
described by N s pointlike particle of mass m, moving in 
the simulation box of size L x x L y x L z . In this scheme, 


the evolution consists of series of streaming and multi¬ 
particle collision steps. In the streaming step, the dy¬ 
namics is evolved by Newton’s equations of motion gov¬ 
erned by forces determined from the total potential en¬ 
ergy V (r Nv , r Ns ) of the system®. In the collision steps, 
which occur at time intervals r, the point-like solvent 
particles are sorted into cubic cells with cell size clq. The 
choice of cell size is such that the mean free path A < a$. 
Multiparticle collisions among the solvent molecules are 
performed independently in each cell, which results in 
the post collision velocity of solvent particle i in cell £ 
being given by v[ = + o)^(vi — V^), where is the 

center-of-mass velocity of particles in cell £ and is a 
rotation matrix. In order to ensure Galilean invariance 
for systems with small A, a random grid shift is applied 
in each direction of the simulation box 35 36 . The method 
described above conserves mass, momentum and energy 
and acco unts fo r the hydrodynamic interactions and fluid 
flow fields^Ml, which are important for the dynamics of 
the active fluid vesicle. 


Membrane and solvent interaction 


Theoretical descriptions of t he vesicle motion arises 
from self-diffusiophoresiJ ^ 14 * 37 *, which in this case is 
based on phoretic motion of the vesicle in external in¬ 
homogeneous chemical field. To model this, the vesicle 
is immersed in a fluid, where N s solvent particles in the 
MPC model is divided in to A and B kind of particles. A 
fraction of membrane vertices are considered to be chem¬ 
ically active where the chemical reaction A -> B takes 
place, producing an inhomogeneous concentration of so¬ 
lute molecules on the scale of the characteristic size of the 
particle, which in turn is responsible for the propulsion 
of vesicle (See fig. fTI) . Such chemically active vertices are 
called catalytic (C) vertices and others are Non-catalytic 
(N) chemically inactive vertices. The A and B solvent 
molecular species interact with all the vesicle vertices 
through repulsive Lennard-Jones potentials, 


Vps(r) = 4 ep 


© 12 -0 


6 1 
-) +T 


, r < r c , (4) 


with the cutoff distance r c = 2 1 / 6 <t. We use the notation 
Vpsi where S = C,N and (3 = A, B to denote various 
interactions between solvent and membrane. In particu¬ 
lar, we take Vac = Vbc — Van ? which are characterized 
by the energy and distance parameters ca and a, respec¬ 
tively; however, interactions between the N vertex and 
B molecules, Vbn, have energy parameter e#. Hence 
the B molecules produced in the reaction on catalytic 
sphere, react differently with catalytic and noncatalytic 
monomers and this asymmetry is an important elem ent 
in the self-propulsion mechanism for the vesicld^H^. In 
order to mimic the fluxes of reactive species into and out 
of the system to drive it out of equilibrium and lead to es¬ 
tablishment of a nonequilibrium stead state, B molecules 
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FIG. 1. A chemically active vesicle composed of catalytic (yel¬ 
low) and non-catalytic (blue) vertices. The figure depicts the 
chemical reaction that leads to conversion of fuel A (orange) 
to product B (green) molecule within the boundary layer, 
hence producing the inhomogeneous distribution of chemical 
species. The boundary layer around the vesicle within which 
the inter molecular forces act is also shown. 


are converted back to A molecules when they diffuse far 
enough away from the vesicle. 


Simulation parameter 

Other simulation details are as follows: all quantities 
are reported in dimensionless LJ units based on energy 
€a, mass m, and distance a parameters. The vesicle move 
in a solvent of A and B molecules within a cubic box of 
volume V and linear dimension 30 with periodic bound¬ 
ary conditions. The MPC rotation angle was fixed at 
4> = 90°. The average number density of solvent out¬ 
side the vesicle was kept at p out ~ 10 , whereas the inside 
solvent number density pi n was varied from 1 to 10 to 
vary shape of the vesicle. The masses of both A and B 
species are m = 1 and the masses of the membrane ver¬ 
tices were adjusted according to their diameters, d v = 2 cr 
to ensure the density matching between vertex and sol¬ 
vent. The system temperature was fixed at T = 0.26. 
Newton’s equation of motion is integrated using velocity 
Verlet algorithm, with MD time step At = 0.005. The 
MPC time interval is r = 0.25 and the time interval for 
bond flip is tbf = 0.1. We have probed various sizes of 
the vesicle where N v varies from 110 to 677. Unless spec¬ 
ified the value of k = 5 and 7 = 80. For the bond and 
excluded volume interaction, we use the cutoff lengths to 
be r c 0 = 1.15 and r c \ = 0.85. The LJ parameters were 
chosen to be ca = 1 . 0 , while e# varies from 0.01 to 0 . 2 . 


III. MEMBRANE VESICLE PROPULSION 

We investigate the self-propulsion of vesicle of differ¬ 
ent shapes- spherical, ellipsoidal and discoidal vesicles. 
These various shapes are generated by changing the num¬ 
ber of solvent particle density enclosed by the vesicle 
which in turn affect the area-to-volume ratio leading to 


(a) 


(b) 


(c) 


9mm 


v = 1.0 


v = 0.85 


v = 0.65 


FIG. 2. Shape transformation of membrane vesicle as a func¬ 
tion of the reduced volume v; (a) sphere (b) ellipsoid and (c) 
discoid 


a shape transformation. The shape characterization is 
done in terms of dimensionless variable, the reduced vol¬ 
ume t>, defined by v = 4^7/3 , where Ro is the radius of 
a sphere with the same surface area as that of the vesicle 
and V ves is the volume of the vesicle. For v = 1, the shape 
is necessarily spherical, which we obtained by keeping in¬ 
side and outside solvent number density same. On reduc¬ 
ing inside solvent density we get the prolate/ellipsoidal 
shape at reduced volume v ~ 0.85, further reduction of 
reduced volume to v ~ 0.65 leads to discocytes vesicles. 
Fig. [ 2 ] shows different shapes obtained in our simulations 
at different reduced volume. 


Spherical vesicle 

The self-propulsion of a membrane vesicle occurs as a 
result of the nonequilibrium concentration gradient pro¬ 
duced by the chemical reaction at the catalytic vertices of 
the vesicle, in conjunction with the different intermolec- 
ular forces between the A and B chemical species and 
the vertices. The ability of the self-propelled vesicle to 
execute directed motion is determined by monitoring the 
mean value of the center of mass velocity of the vesi¬ 
cle motor, projected along the axis between the center 
of mass of the vesicle to the center of mass of catalytic 
vertices: V z = ((V(£) • R (£))), where the double angu¬ 
lar bracket denote the average over time and realization. 
R (t) is the unit vector along the above said axis. An 
isolated vesicle motor undergoes self-propelled motion in 
the given direction with average velocity (V z ); however, 
such small motors are also subjected to strong thermal 
fluctuation that lead to a distribution of propagation ve¬ 
locities. This distribution has been shown to be cl osely 
approximated by a Boltzmann distribution 1 ®!® with 
mean (V z ). The magnitude and nature of the directed 
motion is strongly influenced by the vesicle solvent inter¬ 
actions. Fig.[3]shows the normalized probability distribu¬ 
tion P(y z ) for the system with N v = 302 with Cf = 25% 
(Cf is percentage fraction of total vertex that are cat¬ 
alytic) for different values of energy parameter 

As stated before asymmetry is a crucial component for 
the self-propulsion for the vesicle, in the current system 
mechanism this asymmetry was supplied by the existence 
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FIG. 3. Plot of P{V Z ), the probability density of V z for vesicle 
with N v m 302. Here, the energy parameter ca = 1.0 is used 
for both the curves. Green and pink solid lines are maxwell 
boltzmann fit to data. ( V z ) = 0.0085 for e B — 0.1 and ( V z ) — 
0.0176 for e B = 0.01 


of a strong product molecule B concentration gradient 
and difference in interaction potential of the product B 
molecule with catalytic and non-catalytic vertex. In fig. [ 3 ] 
we calculate the average directed velocity of the vesicle 
(V z ) for interaction potential strength cb equal to 0.01 
and 0.1 keeping ca = 1.0. The figure clearly illustrate 
better propulsion for larger difference of the interaction 
strength i.e. for higher asymmetry. 

Another important element for vesicle propulsion is the 
local steady state concentration of species B in the vicin¬ 
ity of the catalytic C vertices. In order to have a vesicle 
with maximum possible directed movement it is impor¬ 
tant to have a optimum number of catalytic vertices on it. 
Having too many or too less catalytic vertex will not be 
able establish the sufficient nonequilibrium B concentra¬ 
tion gradient around the vesicle; hence, directed vesicle 
propulsion will not be possible. Fig. [4] (Top) shows the 
average velocity of vesicle (V z ) as a function of percent 
fraction of catalytic vertices Cf. From the figure it is 
clear that Cf ~ 25% catalytic gives maximum directed 
velocity (V z ) ~ 0.018 for N v = 302, which are compa¬ 
rable to the previously reported sphere dimer motors 39 . 
This is further supported by higher strength of the local 
steady state concentration gradient of species B in the 
vicinity of the catalytic vertices shown in Fig. [4] (Bot¬ 
tom), for Cf = 25% and 50%. In figure r is the distance 
of solvent B molecules from the center of mass of cat¬ 
alytic vertices. 

The nanoscale objects will not simply move ballisti- 
cally in a given direction due to strong fluctuations from 
surrounding medium 12 . Instead, the motion will be bal¬ 
listic at short times, but at longer times, the motion will 
change to random walk, in which the spans of directed 
motion will be interrupted by random change in direc¬ 
tion. The effective diffusion coefficient can then be cal- 




FIG. 4. Top: Variation of average directed velocity of the 
vesicle ( V z } as a function of percentage of catalytic vertices 
Cf for N v m 302(red circles) and N v = 677 (green square). 
Bottom: The local steady state concentration of B particles 
around the catalytic center of mass ,n B (r), for Cf = 25% and 
50%. 


culated by the velocity autocorrelation function as 

1 r°° 

D e = - d j 0 dt(V(t)-V), (5) 

where V is the center of mass velocity of the vesicle mo¬ 
tor in d-dimension. The center of mass velocity can be 
decomposed in the averaged directed velocity in the R (t) 
direction and fluctuations as, V = R (t)(V z ) + S giving 

D e =D 0 + ^(V z ) 2 T R . (6) 

The first term on the right hand side is the diffusion 
coefficient in absence of any propulsion, while the second 
term is characterized by the decay of the orientational 
correlation function with an orientational relaxation time 
tr. The diffusion coefficient can equivalently determined 
from the mean square displacement (MSD) A L 2 (t) = 
{\ r CM(t) — rcM(0)| 2 ) as D e = dt. There 

are two relevant characteristic time scales in the sys¬ 
tem!^. First is the characteristic diffusion time of the 
product solvent particles around the vesicle given by 
r D = R v /D, where D is the diffusion coefficient of solvent 
molecules and R v is the vesicle size. This time scale sets 
the relaxation time of the redistribution of the particles 
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FIG. 5. Top: Mean square displacement A L 2 (t) for different 
sized vesicles having Cf — 25% depicting the ballistic motion 
~ t 2 at short times and diffusive motion ~ t at long times. 
Bottom: Orientational autocorrelation function for two vesi¬ 
cles. The dashed line is an exponential fit to data. 


around the vesicle when it changes orientation. Second 
is the orientational relaxation time tr which controls the 
changes in the orientation of the spherical vesicle, and 
is defined via the orientation autocorrelation function: 
C e (t) = (R(i) • R(0)> ~ e TR \ The solvent trans¬ 
port coefficient can be calculated analytically for MPC 
dynamicd^, giving tr ~ 600 for vesicle with N v = 110. 

For short times (t <C td) MSD is shown to exhibit 
ballistic regime 40 where A L 2 (t) ~ ( V z ) 2 t 2 , while for 
long times t tr it exhibits a linear behaviour as 
A L 2 (t) ~ 6(D 0 + (V z ) 2 tr/ d)t. Fig. [5|top) shows the 
MSD for the vesicle of different sizes N v = 110 and 302. 
Fitting of the short and long time of this plot is com¬ 
parable with the theoretical expression given above to 
estimate the average vesicle propagation velocity ( V z ) 
and reorientation time tr and enhanced diffusion coef¬ 
ficient. Fig. [^bottom) shows the orientation autocorre¬ 
lation function for different sized vesicle. Table U shows 
the variation of these quantities by varying the size of 
vesicle. 


It is instructive to consider the values of these quanti¬ 
ties for motors with various size. For most of the vesicle 
sizes reported in table [I] we find the diffusion coefficient in 


N v 

(Vz) 

tr 

D e 

29 

0.021 

175 

0.020 

110 

0.019 

1100 

0.132 

194 

0.018 

2500 

0.270 

302 

0.017 

6000 

0.578 

434 

0.016 

11500 

0.981 


TABLE I. The average directed velocity, orientational relax¬ 
ation time and the enhanced diffusion coefficients, obtained 
from the mean square displacement and orientation autocor¬ 
relation function for different sizes of spherical vesicles. 

absence of propulsion to be Dq ~ 10 -3 , whereas the con¬ 
tribution from self-propulsion in D e is much higher. Con¬ 
sequently, the dynamics for these vesicles are dominated 
by the directed motion. However, as decrease the vesicle 
size the directed motion contribution given by \{V z ) 2 tr 
keeps decreasing and, for very small vesicle motors with 
N v = 29, D e is much smaller. This is due to the fact that 
the translational Brownian motion contribution starts to 
dominate and directed motion is compromised. One can 
therefore infer here that the orientational Brownian mo¬ 
tion is not only stochastic element, rather it alters the 
dynamics of chemically powered motor significantly. 

Prolate/Ellipsoidal vesicle 

Vesicles are highly adaptive structures having a rich 
diversity of shapes. Some important reasons that drives 
the sh ape tr ansition are temperature, pressure, curvature 
energy 41 42 -etc. In order to get an insight on the propul¬ 
sion efficiency of vesicle in various shapes, we induce the 
shape change either by changing the number density of 
solvent inside the vesicle, which in turn changes the re¬ 
duced volume v and hence the shape of the vesicl^H or 
by varying the bending elasticity and thus by curvature. 

To evaluate the morphological changes of the vesicle 
in a quantitative manner, we calculated ” asphericity”, a 
as an order parameter that reflects the deviation from 
the spherical shape 31 . Asphericity a is obtained from 
the three eigenvalues (Ai, A2, A3) of the moment-of-inertia 
tensor given as 

(Ai — A 2 ) 2 + (A 2 — A 3 ) 2 + (A 3 — A1) 2 

“=-a?*- <7> 

where, R 2 = Ai+A 2 + A 3 . Fig.pfTop) shows the increase 
in asphericity on reducing the oending rigidity k of the 
vesicle. The observed shape transition in this case is from 
spherical vesicle to prolate/ellipsoidal vesicle. This shape 
change is accompanied with a strong change in dynam¬ 
ical behaviour of the vesicle. Fig. [^Bottom) shows a 
motile prolate which reorients in due course of time due 
to thermal fluctuation. Fig. [7] shows the comparison of 
the probability distribution of directed propulsive veloc¬ 
ity for spherical and ellipsoidal vesicle. It is quite evident 
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FIG. 6. Top: Variation of asphericity a and the average di¬ 
rected velocity ( V z } as a function of bending rigidity k keeping 
pin = 1. Bottom: The propulsion and reorientation of a pro¬ 
late vesicle. 


that the prolate/ellipsoidal vesicle has enhanced propul¬ 
sion velocity, (V z ) = 0.023 as compared to the spherical 
shape, (V z ) = 0.017. It must be noted here the motion 
of the elliptical vesicle is such that the catalytic vertices 
are around the tip of the ellipse and movement is along 
the long axis. This orientation can be explained in terms 
of reduced bending which allows the solvent interactions 
to deform the vesicle thus reducing the drag. 



FIG. 7. Plot of P(\ 4), the probability density of V z for vesicle 
with N v = 302 with prolate and spherical shape. Green and 
pink solid lines are maxwell boltzmann fit to data. 


The increased propulsion velocity in case of pro¬ 
late/ellipsoidal vesicle can be attributed to the less drag 


experienced by the vesicle from its environment in this 
shape. The strength of the drag force depends not only 
on the viscosity at low speeds, but also on the cross- 
sectional shape that is presented to the fluid by the ob¬ 
ject in its direction of motion. For a spherical vesicle of 
radius R the drag coefficient is given by the Stoke’s law 
as ft = 6 ttt] i?, where p is the viscosity of the medium. 
For an ellipsoidal case the drag coefficient changes to 
ft = i n ( 2 a/b)—i /2 f° r m °ti° n at low speed parallel to the 
long axis of the ellipsoid^. Here a and b are semi-major 
and semi-minor axis respectively. The calculated drag 
coefficient for the vesicle having N v = 302 in prolate and 
spherical shapes are found out to be « 98 and 125 re¬ 
spectively. We believe that, it is this decrease in drag 
coefficient which results in enhancement of the directed 
velocity of the prolate vesicle. 

Does the increased directed velocity of the ellipsoidal 
vesicle guarantee it to be a better candidate for mi¬ 
croscale swimmer than a spherical vesicle? To answer 
this next we probed the mean square displacement and 
orientational angular correlation for the prolate shape 
and compared it with the sphere. 




FIG. 8. Top: Orientational autocorrelation function as a 
function of time for ellipsoidal and spherical vesicle having 
N v — 302. Bottom: Mean square displacement A L 2 (t) for 
ellipsoidal vesicles having N v = 302 depicting the ballistic 
motion ~ t 2 at short times and diffusive motion ~ £ at long 
times. 


Fig. [8] (Top) compares the orientation autocorrelation 
function Co(t) for the two different shapes discussed so 
far, which can be used to calculate the orientational re¬ 
laxation time tr. The faster decay of autocorrelation 
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function in case of ellipsoidal vesicle clearly depicts the 
faster reorientation of the vesicle in prolate shape. The 
relaxation time tr for the prolate shape is found out to 
be 1700, which is much smaller than the spherical vesi¬ 
cle having same number of vertices. Table [T] shows the 
spherical vesicle of N v = 302 has tr = 6000. Further, 
the change of directed motion to random walk at the time 
scale of tr is also shown in the fig. [^Bottom). The reduc¬ 
tion in the reorientation time for ellipsoidal case could be 
due increase in floppiness in the prolate vesicle and hence 
that of the catalytic tip. In case the prolate vesicle would 
have been a rigid object, then the concentration gradient 
produced by the vesicle that is responsible for the self 
propulsion would have been along the vector R(t) and 
would have been symmetric normal to the vector. How¬ 
ever, presence of floppiness breaks this symmetry and 
induces forces that lead to reorientation of the vesicle. 


Discocyte vesicle 


As stated before, changing the reduced volume v leads 
to various shapes of the vesicle. Next we investigate the 
propulsion of vesicles having a discocyte shape, which 
was obtained by reducing the number density of solvent 
inside the vesicle to pi n = 2 with k = 5 . Unlike the case 
of ellipsoids, where we have a unique mode of propulsion 
with the velocity direction along its long axis, here we 
obtain two different modes of propulsion, in one the dis¬ 
cocyte moves with its face on, and in the other it moves 
with its edge at the front. Fig. fright panel) shows both 
the propulsion configuration. Fig. [9] (left panel) compare 



FIG. 9. Left Panel: Comparision of P(V z ), the probabil¬ 
ity density of V z for vesicle with N v = 677 in discocyte and 
spherical shape. Right Panel: Instantaneous configuration of 
moving discocytes in the face on and edge-on mode. Arrow 
shows the direction of movement of the vesicle. 


the directed velocity V z for the discocyte with a spher¬ 
ical vesicle having same vertex number N v = 677. It 
is evident from the plot that the discocyte have higher 
propulsion with (V z ) = 0.019 than its spherical counter¬ 
part with (V z ) = 0.014. 


IV. VESICLE MOTOR EFFICIENCY 

Chemically powered vesicle motors converts chem¬ 
ical energy into mechanical work, driving the self- 
propelled directed motion, while working in an environ¬ 
ment where viscous drag is important. The efficiency r] s 
of chemically-powered motor, which measures the effec¬ 
tiveness of mechanochemical energy transduction, can be 
defined a d 11 * 45 * 


Pmech . C(V) 2 

Pzhem AfJ,R 


( 8 ) 


where, P m ech is the mechanical power output of the 
motor and P c hem is the total chemical power input to 
the motor. Pmech and P c hem are further defined as 
Pmech = C(V) 2 and Pchem = AfxR, where £ is the 
drag coefficient of the vesicle motor, (V z ) is the directed 
propulsion velocity of the vesicle, R is the net chemical 
reaction rate and A/jl is the change in the chemical po¬ 
tential in the reaction. To determine the efficiency from 
above equation, the net chemical reaction rate R was 
calculated by counting the number of A B reactive 
events that take place at the catalytic vertices per unit 
time. In the MPC-MD simulations the change in the 
chemical potential is given by 


A/i = p B ~ Pa = -K B T In 


n B nl 

eq 5 

n B n A 


(9) 


where n ^ and n e ^ denote the equilibrium number den¬ 
sities of A and B species respectively, while n A and n B 
are the steady state densities^. The drag coefficient de¬ 
pends not only on the viscosity, but also on the cross- 
sectional shape that is presented to the fluid by the vesi¬ 
cle in its direction of motion. For spheres, ( = 67 rryR, 
for ellipsoids moving along its long axis ( = 1/2 

and for disks moving edge on £ = ^ 77 p, where p is the 
radius of the disk. We compare the Stokes efficiency of 
the various sizes and shapes of the vesicle in Table [H 


N v 

Shape 

e# 

<v> 

Vs 

110 

Spherical 

0.01 

0.019 

0.0027 

302 

Spherical 

0.01 

0.017 

0.0020 

302 

Spherical 

0.1 

0.008 

0.0004 

302 

Ellipsoidal 

0.01 

0.023 

0.0030 

677 

Discocyte 

0.01 

0.019 

0.0015 

677 

Spherical 

0.01 

0.014 

0.0012 


TABLE II. Comparison of the average directed velocity and 
Stoke’s efficiency for vesicle motor of various size and shapes. 
The LJ energy parameter e A = 1.0. 


It is evident from the above calculation that the power 
associated with the mechanical work of the motor is 

















larger for smaller vesicle of same shape, as the Stokes 
efficiency for spherical vesicle of N v = 110 is 0.27% and 
that of vesicle with N v = 302 is 0.20%. On changing the 
shape of the vesicle its efficiency also changes. Table [Tl 
shows that the ellipsoidal motor with N v = 302 is more 
efficient with rj s = 0.30% than its spherical counterpart 
with T] s = 0.20%. Similarly for comparison of spherical 
and Discocyte shape shows enhancement of efficiency in 
the Discocyte shape. The maximum efficiency achieved 
in our simulations is about 0.30%, which is similar to 
the other chemically powered motonJ^HS however the 
this efficiency is much less than that of most nano- and 
micromotors in biolog}®. 


V. CONCLUSION 

The particle-based mesoscopic model for chemically 
powered self-propelled vesicle dynamics discussed in this 
article provides insight into the factors that control the 
motion. Our model couples the vesicle motion to the sol¬ 
vent by intermolecular forces, which ultimately manifest 
themselves in properties such as the surface tension of 
macroscopic scales. Here, we show that depending on 
the bending rigidity and solvent density around the vesi¬ 
cle, different shapes of vesicle can be obtained and all 
these shapes exhibit very different dynamical behavior. 

Many potential applications of such vesicle motors will 
involve launching them to perform a given task, such as 
cargo transport, and further investigations of interaction 
between them can provide the information needed to de¬ 
sign motors that can cooperate with each other to per¬ 
form such tasks. 
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